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Abstract: We study scalar excitations of high spin operators in A/" = 4 super Yang-Mills 
theory, which are dual to solitons propagating on a long folded string in AdS^ x S^. In the spin 
chain description of the gauge theory, these are associated to holes in the magnon distribution 
in the s[(2, M) sector. We compute the all-loop hole S-matrix from the asymptotic Bethe ansatz, 
and expand in leading orders at weak and strong coupling. The worldsheet S-matrix of solitonic 
excitations on the GKP string is calculated using semiclassical quantization. We find an exact 
agreement between the gauge theory and string theory results. 



1 Introduction 



The discovery of integrability on both sides of the AdS/CFT correspondence has opened the 
possibihty for precise tests of the conjecture. On the gauge theory side, the full spectrum of 
planar A/" = 4 super Yang-Mills theory has been determined by the all-loop asymptotic Bethe 
ansatz [1] for an integrable p5u(2,2|4) spin chain. Of particular interest is the si{2) subsector 
of twist operators composed of light-cone covariant derivatives acting on complex scalar 
fields Z, 

tr(ZD^Z) + • • • , (1.1) 

where the ellipsis denotes all other operators that mix under renormalization. At large Lorentz 
spin 5*, the conformal dimension A of the ground state in this sector exhibits logarithmic scaling 
at all loops 

E = A-S = 2reusp(^7) log S + 0(log° S). (1.2) 

The cusp anomalous dimension Fcusp interpolates smoothly between small and large values 
of the coupling constant g = VX/Att [2]} It also controls the infrared divergences of gluon 
scattering amplitudes, which are dual to light-like polygonal Wilson loops. There has been 
much recent interest in studying excitations around this ground state [4-6] . In particular, their 
spectrum is an important ingredient in the operator product expansion of light-like Wilson 
loops [7]. The spectrum of all such excitations has been determined by Basso [8], by way of 
a set of integral equations for the excitation densities which can be solved at all values of the 
coupling. The scalar excitations on top of the sea of derivatives (1.1) are known as holes. In 
the spin chain language, the holes are gaps in the magnon distribution. The vacuum contains 
two large holes of order 0{S) which contribute to the logarithmic scaling [9]. Excitation of the 
vacuum introduces additional holes of order 0{S^). 

The string state dual to the ground state operator of fixed spin (1.1) is a long folded 
string rotating in AdS^ x 5*^, known as the GKP string [10]. The string energy also exhibits 
logarithmic scaling due to the length of the string in the target space, whose coefficient coincides 
with Fcusp at strong coupling. The two spikes at the end of the string approach the boundary 
of AdS^ and are identified with the two large holes. More generally hole excitations correspond 
to classical solitons that propagate on the GKP string. These are identified with solitons of the 
sinh-Gordon equation after Pohlmeyer reduction of the string equations of motion [11]. The 
classical dispersion relation for these excitations was obtained in [12] where they were named 
"Giant Holes" in analogy with the Giant Magnons of [13]. 

In this paper we study the scattering of two holes/solitons for all values of the coupling. 
On the gauge theory side, the scattering phase is obtained by considering the scattering of 
the second hole with the change in magnon density due to the first hole [14]. Following the 
techniques introduced in [8, 15], we obtain an integral expression, which we expand at weak and 
strong coupling. On the string theory side, the worldsheet S-matrix is calculated from the time- 
delay [16] that one soliton accumulates as it passes the other. It receives additional contributions 

"'^For a review of twist operators and the cusp anomalous dimension, see [3]. 
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due to the non-trivial relation between worldsheet coordinates and global coordinates and also 
due to the contribution of each excitation to the length of the string. Once these subtleties are 
carefully taken into account, we find exact agreement with the gauge theory result at strong 
coupling. 



2 Scattering of spin chain holes 

First it is important to be precise about we mean by hole excitations. Consider operators of 
arbitrary twist 

tr(D^Z^) + ■ ■ • . (2.1) 

In the spin chain description, we identify the S covariant derivatives with magnons propagating 
on a background of scalars corresponding to the ferromagnetic ground state. An important 
feature of the non-compact s[(2,M) spin chain is that there can be an arbitrary number of 
magnons at each spin site without increasing the length of the chain which is identified with 
the number J of scalars. In addition, the number of holes in the mode number distribution of 
the magnons is also equal to J. The Bethe ansatz equation provides a quantization condition 
for the magnon rapidities. The large spin limit we will consider represents a highly excited state 
in which the total length is held fixed. In this limit, it is natural to view the sea of derivatives as 
the pseudovacuum on which the scalars act as hole excitations. Moreover, Alday and Maldacena 
have identified the elementary fluctuations around the GKP string with insertion of twist-one 
operators in the pseudovacuum [4]. With these motivations, we define the holes as vacancies 
in the magnon distribution, to be identified with inserting a scalar Z, as opposed to removing 
a derivative in the trace (2.1). 

2.1 One-loop S-matrix at weak coupling 

The one-loop Bethe ansatz equation for 5* magnon excitations in an s[(2,]R) spin chain of J 
sites is 

s 

l = e-^P^^'^^'l[S{uk,ue), (2.2) 

where the magnon momentum and S-matrix are defined as 

p{uk) = ilog — — , S{uk,U(,) = — . (2.3) 

Uk + ij 2, Uk — U£ + I 

Each Bethe root Uk is associated with a mode number n{uk) corresponding to a branch of the 
logarithm 



s 



= 27rm Mfc -2 > log — - Jlog — — . 2.4 

Uk-Ui + i Uk + 1/2 

The roots all lie on the real axis. When J = 2, the mode numbers are symmetrically distributed 
as n{uk) = k, k = ±| . . . , ±^^^.^ In the large spin limit, the magnon roots form a cut in 



^For arbitrary J, the magnon mode number develops a gap in the middle corresponding to J — 2 small holes. 
The gap closes in the large spin limit. 
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the interval [— f , f ] and (2.2) becomes an integral equation for the magnon density p{u) = ^ 
[17] 

/2 o(v) 4 
dv ^ -. (2.5) 
_| l + (M-t;)2 l + 4u2 ^ ^ 

Introducing an excited hole of rapidity Uh shifts J = 2 i— )■ 3 and the magnon mode number 
n{uk) I—)- n{uk) + 0{u — Uh), where 9{u) is the Heaviside step function. The change in the ground 
state density po{u) due to the hole, ph{u,Uh), can be written as 

Po{u) + Phiu, Uh) = p{u) + S{u-Uh). (2.6) 

We will also use Ph{u) when we do not need to specify the hole rapidity. It satisfies the integral 
equation 

J 1 + [u — vj^ 1 + [u — Uh) 1 + 4u"^ 

Note that we have extended the integration domain to the entire real line. This is justified 
because the excitations are of order C(S'°) and Ph{u) is suppressed at large rapidity. We solve 
for Ph{u) by Fourier transform to obtain 

Ph{u) = [Ir + iu] (l: -iu] - {1 + i{u - Uh)) - ^ (1 - «(« - Uh)) , (2.8) 



2 J \2 

which decays as ph{u) = ^ + C(m^^) at infinity. 

The S-matrix between two holes carrying rapidities Mh,i,M/i,2 is a phase S{uh^i,Uh,2) = 
— exp(i(5spin). It is given by integrating the scattering of the second hole with the magnon 
density change introduced by the first hole [14]. 

5„^2 \ogS{uh,i,Uh,2) = J du dulogS{u,Uh,2) {ph{u,UhA) - 5{u - Uh,i)) . (2.9) 

The integral can be evaluated by Fourier transform. We find 

T{1 - i{uh,i - Uh,2)) ., T{l/2 + iuh,i) .. T{l/2-iuh,2) 

Sspin = - « log , - t log :^^y:r77^ — ^ - dog 



T{l + t{uh,i~Uh,2)) ^T{l/2-tUh,i) ^T{l/2 + tUh,2) (2.10) 

-P{uh,l) +p{Uh,2), 

where in recovering the phase 5spin from its derivative (2.9) we fixed the constant of integration 
such that the expression is antisymmetric in Uh,i and Uh^2, as is necessary for a unitary S-matrix. 
Note the appearance of momentum terms p{uh,i) and p{uh,2)- They are present because by 
introducing a hole we have created a vacancy in the magnon distribution. In defining the 
effective Bethe ansatz for holes, it is more appropriate to absorb such terms as part of the 
phase due to propagation instead of as part of the S-matrix. ^ Assuming factorized scattering. 



^We thank Benjamin Basso for insightful comments on interpreting the S-matrix and its connection to the 
effective Bethe ansatz. 



- 4 - 



we may write down the effective Bethe ansatz for a twist J chain 

J 

l = \{S{uh,i,Uh,j). (2.11) 

The two large holes Uh,i, Uh,j of order 0{S) are non-dynamical because their positions are fixed 
by the quantum number of the spin chain. We separate their contribution from the other 
dynamical holes. Using Stirling's approximation, we find 

1 = g-.p.K,.)x21ogS+O{50) -Q S{uH,,Uh,,), (2.12) 

i=2 

where Ph{uh) = '^Uh + 0{g'^) is the hole momentum [9] and we redefine the S-matrix as 

^/ . _ r(i - ijuh,, - Uh,j)) r(i/2 + iuh,^ r(i/2 - iuhj) , . 

^ r(l + z{un, - un,)) r(l/2 - ^u^,) V(l/2 + lun,) ' ^ ^ ^ 

The two large holes define the effective length of the chain (=2 log 5*) as seen by the other 
dynamical holes. There are 0{S^) contribution to the effective length from the momentum 
terms in (2.10) that was absorbed into the propagation phase, as well as an additional phase 

—2i log — — -. — ' from separating the two large holes. They are suppressed at large 5*. The 
T{l/2-iUh,i) 

effective Bethe ansatz (2.12) is equivalent to the quantization conditions of [9] (see equation 
(3.41) of this reference). 

By comparing (2.8) and (2.13), we note that the excitation density ph{u,Uh) coincides with 
the hole scattering kernel idu,^ log S{u, Uh). This is not surprising because the definition of Ph{u) 
(2.6) can be rewritten as an integral equation for the hole density 5{u — Uh) 

Pq{u) - p{u) = 5{u - Uh) - j dv ph{u,v)5{v - Uh). (2.14) 

Hence ph{u) can be naturally interpreted as the hole scattering kernel. This general fact has 
also been observed in the study of Destri-de Vega type non-hnear integral equations [18, 19]. 

2.2 All-loop asymptotic Bethe ansatz 

The all-loop generalization of (2.2) is given by the asymptotic Bethe ansatz equation [1] 

1= ^ n^^"'^'"^)^'''^"""'^' (2.15) 
where the magnon S-matrix is defined by the deformed rapidity -Ujt ± | = + ^ as 



^k -xj K'^- g'^/ixk^i) 



s{u„ u,) = ^-^^ ( ; , (2.16) 



- 5 - 



and 6{u,v) is the BES dressing phase [20] that gives the correct answer beyond three loops in 
weak-couphng and at leading order in strong coupling. As in the one-loop case, we obtain an 
integral equation for the excitation density ph{u) in the large spin limit 

= 2iTph{u) + dv [idulogS{u,v) - 2du0{u,v)] ph{v) 

+ 2du9{u, Uh) - idu log ^(u, Uh) + id^ log I ^ . 

J 

Following [8], we split the density into parity even and odd parts in u as phiu) = —ah{u) — 
Uhiu). We also split the magnon scattering kernel idu\ogS{u,v) into IC°,IC°, the dressing 
contributions into Xdressmg,2dressing, and the inhomogeneous term idu\ogS{u,Uh) into Ih,^h- 
From their integral expressions [8], Basso observed that the main scattering kernels are also 
even and odd in v. Thus the parity even and odd parts decouple and we obtain an integral 
equation for each 



= 2nah{u) + j dv lC{u,v) ah{v) + Xdressing («) +1h{u) 
= 2Tcah{u) + j dv jC°{u,v)ah{v) +idressmg{u) +ih{u) 



(2.18) 



We proceed to solve for the densities and compute the S-matrix between two excited holes. As 
before, the scattering phase is given by integrating the excitation density against the inhomo- 
geneous term 



(2.19) 



du^i^spin = J du Ih{u,Uh,2) (Th{u,Uh,l) +Ih{u,Uh,2) '^h{u,Uh,l 
- ^h{Uh,l, Uh,2) - ^h{Uh,l,Uh,2), 

where the even and odd inhomogeneous terms Th,Xh admit the integral representations [8 

POO 

Ih{u,Uh) = -2 / dt cos{ut)e~^^^'Kh{t,Uh), 
Jo 

^ /"OO 

Ihiu,Uh) = -2 / dt sin(Mt)e"*/^7rft(t,M/,) 



' 5 
(•CO 



.OO (2.20) 

Tihit, Uh) = cos{uht)e-^/^ - M2gt) - Ag^ / ds tK{2gt, 2gs) cos{uhs)e-''^ , 

Jo 

POO 

TTh{t, Uh) = sin(M/jt)e"*/^ - Ag'^ / ds tK{2gt, 2gs) sm{uhs)e~''^'^ . 

Jo 

The symmetric BES kernel K{t, s) can be split into parity even and odd parts K = + 
that admit expansion over Bessel functions as 

2 oo 2 °° 

K+it,s) = -5^(2n-l)J2n-i(t)^2n-i(s), K_{t, s) = - ^ (2n) (t) (s) . (2.21) 

n=l n=l 
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Define the Fourier-Laplace transformed densities as 

Qf^^t) = e"*/^ / du ah{u) cos{ut), Qh{t) = e"*/^ / du ah{u) sm{ut), (2.22) 



for t > 0. By using the integral representations for the main scattering kernel, the dressing 
contributions, and the inhomogeneous terms [8], one can write the integral equations (2.18) as 
equations for Qh, 



(e* - l)Qh{t) + V ^ ds t (^K{2gt, 2gs)Qh{s) + K.{2gt, 2gsf-^^ ^ 



- 1 J 



(2-23) 



(e*-l)fi,(t)+V^ dst(^K{2gt,2gs)nt,{s) + K+{2gt,2gs)'^-^^^ =Mt), 

where , are the generating functions for the conserved charges sourced by the holes. Their 
dependence on the excitation densities can be seen from the following integral representations 
[8] 



POO 

qti2gt) = 8g^ / ds tK+{2gt, 2gs) (^^(s) + cos(ufes)e"'/^) , 
Jo 

ql{2gt) = 8g^ ds tK_{2gt, 2gs) (fi^s) + sin(M;,s)e-^^ 



(2.24) 



One recognizes (2.23) as a generalization of the celebrated BES equation [20], which gives an all- 
loop expression for the density that can be expanded at weak and strong limits of the coupling 
and is amenable to numerical study at intermediate values of the coupling. The expression for 
the scattering phase (2.19) can be written in terms of fi/^, VL^ as 



dufi2^spin 2 / dt 



{Qh(t, Uh,i) + e cos(Mh,it)) 7r/i(t, Uh,2) 

(2.25) 

+ \SLh{t, Uh,i) + e~*/^ sin(u,,,it) j -Khit, Uh,2) ■ 

At weak coupling, we can solve (2.23) iteratively by expanding the Bessel function near the 
origin. We recover the one-loop scattering phase (2.10) from (2.25). 

2.3 S-matrix at strong coupling 

To expand (2.25) at strong coupling, it is useful to rescale u = 2gu and introduce the auxiliary 
densities 7^(t),7^(t) related to ^lh{t),D.fi{t) by 

(e*/^^ - m,{t) = cos{u,t)e-'/'^ - Mt) + ^t(t) + ^X{tl 
{e"^^ - l)fi,(t) = sin(w,t)e-*/^^ + ^^{t) + 7^(t). 

The auxiliary functions 7±,7± admit expansion as Neumann series of Bessel functions 

oo oo 

7^^ = 2^(2n- l)72Vi^2n-i(t), = 2Y,{2nhU2n{t), (2.27) 

n=l n=l 
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and similarly for j!^. We expand (2.23) in series of Bessel functions and separate into parity 
even and odd parts. We obtain four equations, two of which relate the higher conserved charges 
to the auxiliary functions, q'^ = —^■y^, = —2^^, which can be used to express the other 
pair solely in terms of the auxiliary functions 



7n + 



t 

dt 



Jn{t) 



gt/29 _ 1 

l1{t) + (-l)"7^(t) 

pt/2g _ 1 



dt Jn{t) 
t e*/29 - 



(2.28) 



sm{uht)e 



We follow the technique in [15] to transform the set of equations (2.28) into a singular integral 
equation. Define V{t) = T'l{t) + iVlit) related to -f''{t) = 7^(t) + i-f'^{t) via 



r^(t) = 7i(t) T coth(t/%)7i(t), r^(t) = ±%{t) + coth(t/%)7i(t). 



(2.29) 



At strong coupling, the first term in the LHS of (2.28) are suppressed. We apply the Jacobi- 
Anger expansion to replace the infinite system of equations (2.28) with a single equation for 



and 



/oo i-oo iut 

dt [e'^*r^(t) - e-'^'r1{t)] ^ -2 dt ^^75^ [cos(M,t)e*/^^ - Mt)] 



dt 



e*"T!l(t) -e-*"T'!(t) 



(2.30) 



dt 



.t/2g 



sm{uht)e 



with |n| < 1. In the Fourier space, they become singular integral equations for r^(u) and T'^{u) 
^ ' +7rr^(M) ^ -2g [log {A\u^ -ul\) +2 sm-\u)] = -2/€^(m), 

(2.31) 



dk 



/ 



dk 



k — u 



7rsgn(M/,) + log 



U + Uh 



U-Uh 



-2k^{u). 



The above relations hold for \u\ < 1. For \u\ > 1, T^{u),T^{u) can be determined by analyticity 
and is subleading at strong coupling [15]. Thus we can restrict the range of integration to the 
unit interval. Such singular integral equations have a standard solution 



f''(A;) 



A:2<1 1 



TT 



K\k) - 



1 + k 
1-k 



1/4 



^ du k'^' 



u — k 



1-u 
1 + u 



1/4 



V2c 
1 + k 



(2.32) 



where c can be determined from the so-called quantization conditions [15] and is negligible. 
The finite Hilbert transform in (2.32) can be evaluated using contour integration [21]. As 

k'^{u),k'^{u) have branch points, it is simpler to consider instead du,^T'^^{u) and duf^T'^{u). We 
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find 



V2g fl + k 



TT 
TT 



k 



1/4 



1/4 



Uh 



Uh + 1 



1/4 



1/4 



Uh 



Uh + lJ k-Uh 



+ 



Uh - 1 

Uh + l 
Uh-l 



1/4 



k + Uh 



1/4 



k + Uh 



(2.33) 



We may now compute the scattering phase (2.25) from the transformed densities r^{t),r^{t). 
Performing the integral is straightforward but tedious. Remarkably, the answer takes a simple 
form when expressed in terms of the Zhukovsky variables Vi = Uh^i — \j^h 
to the soliton velocities in the string worldsheet 



^ — 1, which correspond 



^spin 



2^? 



(P1E2 - P2^l) + 2(? 



V2 



1 

V2 



\ogv + 





(i 




7i72 




V2J _ 



(2.34) 



where Pi,Ei are the energy and momentum of the giant holes (3.10), v is the soliton velocity 
in the center of mass frame (3.14), and 7^ = (1 — vf)~^^'^ is the Lorentz factor. The phase is 
antisymmetric in Vi and V2, as is expected from a unitary S-matrix. 



3 Worldsheet scattering of giant holes 

In this section we review the construction of solitonic excitations on the GKP string [12] and 
study their scattering. Consider classical strings in AdS^, in embedding coordinates 

-Xo'-X2 + x2 + X| = -l. (3.1) 

We may form two complex coordinates and express in global AdS^ coordinates {p, 0, t} as 

= Xo + iXi = e** cosh p, Z2 = X3 + iX^ = e''^ sinh p. (3.2) 

Integrability of the classical string is manifest by performing a Pohlmeyer reduction on the 
equation of motion to obtain the sinh-Gordon equation [11] 

d^a + sinha = 0. (3.3) 

The GKP string is a vacuum solution with sinh-Gordon angle a = 0, which yields 

= e'^ cosh a, = e'^ sinh a. (3.4) 

As in [12], to regulate the IR divergence in the string energy a cut-off ±A on the range of the 
worldsheet coordinate a is introduced. It is natural to consider solitonic excitations of this 
vacuum. However, a one-soliton solution does not correspond to a physical state of the string 
due to the zero total momentum constraint Pj = 0. Instead, we construct a two-soliton 
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solution and then identify one of the sohtons when it is well-separated from the other as the 
"one-soliton solution". The two-soliton solution in the center of mass frame is [11] 

^ssr IT cosh T cosh cr + cosh X cosh cr — a/1 — f ^ sinh X sinh a — — f ^ sinh T cosh a 

^ ' cosh T + V cosh X ' 

,v V cosh T sinh a + cosh X sinh a — \/l — v"^ sinh X cosh cr — i-\/l — sinh T sinh a 



cosh T + V cosh X 



where X = 2'~fa,T = 2v'JT with sinh-Gordon angle 



a. 



log 



V cosh X — cosh T 



(3.5) 



(3.6) 



f coshX + coshT^ 

We may study a single soliton by considering a limit of the above solution where the second 
soliton is sent to infinity. The solution is obtained by taking the large T, X limit of (3.5) and 
takes the form [12] 

j^e^(cosh(T — a/I — f 2 sinho") + e~^{v — — v"^) cosh a 
j^e^(sinha — a/1 — f ^ cosh a) + e~^{v — za/1 — t>2) sinh a 



r, cr 



^2(r,a) 



where S = (X — T)/2 = 7(0" — fr). The sinh-Gordon angle reduces to that of a single soliton 
so it can be interpreted as a one-soliton solution. At large distance the soliton solution can be 
matched onto the vacuum or GKP solution using the asymptotics, 



as (J — )■ cxo, 



as a — )■ —CO. 



(3.8) 



In order to match the asymptotics of the GKP string solution, we need to restrict the worldsheet 
coordinate to —A < cr < A — a(f ) where A is the IR cut-off [12]. Note also that global AdS 
time no longer coincides with the worldsheet time and is shifted by a constant t ~ r — /3(w) as 
cr —00 where we define, 



a[v] 



log 



1-v^ 



tan 



The dispersion relation of the giant holes has been computed in [12] as, 



E{v) = g 



log 



7 + 1 
7^ 



P{v) = 2g 



7f 



tan 



-1 



1 

7f 



(3.9) 



(3.10) 



To move to a frame where the two solitons have generic velocities vi,V2, we introduce the 
following boost 



cosh^ -sinh^ 
-sinh^ cosh^ 



(3.11) 
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where 6i is the sohton rapidity. One may easily check that 



]{X - T) = 71 (a' - v,t'), bx + T)= ^2{cr' - v^r'). 



(3.12) 



2' / -V 2 

As r' — > oo, the two sohtons are now located at a' = Vir'. We may separate them by zooming 
into the first soliton in the limit T — )■ ±oo and X — )■ ±oo (for the second soliton we would 
zoom in the limit T — t- ±oo and X — t- =Foo). The sinh-Gordon angle will be that of a single 
soliton with a time-shift. 



ttss logtanh^ <^ 7i 



as t' — )■ ±oo. 



(3.13) 



We may read off the time delay experienced by the first soliton from (3.13) as 

logv 2 7172(1-^^1^2)-! 



12 



7ifi 



7172(1 - V1V2) + 1' 



(3.14) 



The scattering phase exp (iS) can be calculated from the semiclassical formula = AT12 [16] , 
where Ei = E{vi) is the soliton energy evaluated at vi. Although E{v) was computed in the 
center of mass frame, it is intrinsic to the single soliton so is independent of frame. We integrate 
the time delay to find 
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72V2 



log 



7i + 1 
7i - 1 



tan 



72 



1 



71^1 



+ 2 Ui 



1 



V2 



1 

V2 



logv 



(3.15) 



where we used the asymptotic condition 5 — )■ as tii — > 1 to fix the integration constant. 
However, the S-matrix does not appear unitary: 6{vi,V2) 7^ —6{v2,Vi). This is familiar from 
the giant magnon case [13] (see equation (3.33) of this reference), where a similar problem arose 
and was related to an ambiguity in the length assigned to the excitation. In our case, the second 
soliton changes the string length by —02 = —C({v2)- Because the energy E = A— S" is canonically 
conjugate to (t + 0)/2, we also need to include the additional time delay (32 = /3(f2) from the 
second soliton due to the difference between the global AdS time t and the worldsheet time r. 
Taking these effects into account, the scattering phase is given as ^string = S — Pia2 + Eif32. 
It is then antisymmetric in Vi and V2 and exactly matches the gauge theory result at strong 

coupling, i.e. ^string = (^spin- 



4 Discussions 



In this paper we have considered the scattering of two holes on a spin chain of length two. The 
generalization to arbitrary number of holes on a chain of arbitrary length is straightforward. 
The string theory picture is that of multiple solitons propagating on a folded, more generally, 
spiky string [22]. Classically we know that these solitons undergo factorized scattering and, as 
the theory is integrable, it is natural to expect that the S-matrix for any number of holes is 
exactly factorizable for all values of the coupling. In principle, the two hole S-matrix obtained 
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in this paper, together with the dispersion relation of [12], can therefore be used to write down 
a dual set of asymptotic Bethe ansatz equations which determine the full spectrum of excited 
states of a GKP string of large but finite length (in the non-compact SL{2) sector). This must 
be, of course, equivalent to the original all-loop asymptotic Bethe ansatz equations of [1] but 
should be a simpler starting point for discussing generic states of large spin. 
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